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ABSTRACT 

We present observations of the BL Lac object 1ES 0414+009 in the >200 GeV 
gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was 
observed by VERITAS between January 2008 and February 2011, resulting in 56.2 
hours of good quality pointed observations. These observations resulted in a detection 
of 822 events from the source corresponding to a statistical significance of 6.4 standard 
deviations (6.4 a) above the background. The source flux, showing no evidence for 
variability, is measured as (5.2±l.l s t a £ ± 2.6 sys )x 10~ 12 photons cm -2 s _1 above 200 
GeV , equivalent to approximately 2% of the Crab Nebula flux above this energy. The 
differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with an 
photon index of T= 3.4 ±0.5 s t a t±0.3 sys and a flux normalization of (1.6±0.3 S f a i±0.8 S3/s ) 
x 10~ n photons cm -2 s^ 1 at 300 GeV. We also present multiwavelength results taken 
in the optical (MDM), X-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these 
results to construct a broadband spectral energy distribution (SED). Modeling of this 
SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe 
emission from the system, with a lepto-hadronic model providing a better fit to the data. 

Subject headings: BL Lacertae objects: individual (1ES 0414 +009, VER J0416+011), 
Gamma rays: galaxies 



1. Introduction 

Among the more substantial accomplishments of the latest generation of imaging atmospheric 
Cherenkov telescopes (IACTs) is the reliable detection of increasingly distant active galactic nuclei 
(AGN). While the previous generation of ground based gamma-ray experiments such as the Whipple 
10 meter telescope and HEGRA detected comparatively close AGN, the furthest being H 1426+428 



at a redshift of 0.129 (jHoran et al.1 12002). observations with VERITAS, HESS, and MAGIC have 



27 shown that it is possible to det ect sources of very h igh energy (VHE) gamma rays (>100 GeV) as 



far out as redshifts of z= 0.54 (jAlbert et al.ll2008al ). 



29 The AGN 1ES 0414+009 is among the furthest detected TeV BL Lacertae objects, having a well 

30 measur ed redshift of z=0.287 (jllalpern et al.lll99lh . 1ES 0414+009 w as first discovered in an X-ray 



31 survey performed by the HEAO 1 A-l instrument (jUlmer et al.lll980l ) and was originally associated 



32 with a cluster of galaxies. Further radio, X-ray and optical observations led to the association 



19716, USA 

29 Physics Department, Columbia University, New York, NY 10027, USA 
30 Argonne National Laboratory, 9700 S. Cass Avenue, Argonne, IL 60439, USA 
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of the X-ray sour ce with an X-ray bright BL Lacertae obje ct (jUlmer et al.lll983l ). Additional X- 
ray observations (|Wolter et al.l Il998l ; iBeckmann et al.1 12002 ) show 1ES 0414+009 to have typical 
spectral characteristics of "high frequency peaked BL Lacertae" (HBL) objects, as well as being 
comparatively bright in the 2-10 keV band with a luminosity comparable to Markarian 421 and 
PKS 2155-304, allowing the possibility that 1ES 0414+009 may be a high-redshift analog of these 
well studied TeV blazars. 

1ES 0414+009 also appears in the Fermi-LAT 2-year source catalog, associated with the source 
2FGL J0416.8+01o£] Its inclusion in the 2FGL catalog is based on its detection in the 1-100 GeV 
band at a significance of 6.8 a; the photon spectrum of the excess is well fit by a power law with a 
photon index of 1.96+0. lQ s tat and an integral flux in the 1-100 GeV band of (6.9+1.4 s t a j )xl0" iu 
photons cm -2 s _1 . 

1ES 041 4+009 was among the 33 objects deemed to be very good candidates for TeV detections 
in the list of ICostamante fc Ghisellinil ([2002]) and, as such , preliminary observations were taken by 
several TeV observa tories with only upper limits resulting ([Perez et al.ll2003l ; lAharonian et al.ll2004l ; 
Albert et al.ll2008bl ) . The HESS array of IACT telescopes detected the source in observations taken 
from 2005 to 2009 ^Hoffman et al.ll2009h . with 224 excess events above 200 GeV in 73.7 hours o f 
observation, corresponding to a statistically significant excess of 7.8c (jAbramowski et al.ll2012l ). 



The HESS excess is well fit by a power law with a photon index of 3.5±0.3 s t a t + 0.2 SJ/S and an 



integral flux above 200 GeV of (1.9±0.2 stat ± 0A sys ) x 10" 
Crab Nebula flux in the same range. 



-12 



photons cm 



- 2 s- 1 or ~0.6% of the 



VERITAS Observations 



54 The VERITAS array (jOng et al.ll2009l ) of IACTs located in southern Arizona (1.3 km m.a.s.L, 

55 31°40'30"N, 110°57'07" W) began 4-telescope array observations in September 2007 and is the most 

56 sensitive IACT for observations above 200 GeV currently in operation. The array is composed of 

57 four 12m diameter telescopes, each with a Davies-Cotton tessellated mirror structure of 345 12m 

58 focal length hexagonal mirror facets (total mirror area of 110 m 2 ). Each telescope focuses Cherenkov 

59 light from particle showers onto its 499-pixel PMT (photomultiplier tube) camera. Each pixel has 

60 a field of view of 0.15°, resulting in a camera field of view of 3.5°. VERITAS has the capability to 

61 detect and measure gamma rays in the 100 GeV to 30 TeV energy regime with an energy resolution 

62 of 15-20% and an angular resolution of <0.1° on an event by event basis. 

63 Between June and September 2009, VERITAS underwent a significant reconfiguration in which 

64 one of the four telescopes was moved to a new location making the geometr ic configuration of the 

65 array more symmetric and increasing the overall baseline (jPerkins et al . 2009). This r econfiguration 



66 as well as a marked improvement in the process for VERITAS mirror alignment (jMcCann et al 
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67 12010 ) . yielded a significant improvement in the VERITAS sensitivity at the level of 30%. In its 

6B current configuration, VERITAS can detect a 1% Crab Nebula flux at a 5<r significance in under 30 

69 hours, improved from the ~45 hours required with the previous array configuration and alignment. 

70 The VERITAS observations of 1ES 0414+009 were made between January 2008 and February 

71 2011 and therefore comprise data taken in both array configurations. After quality selection cuts 

72 which remove data taken during poor weather and hardware conditions, approximately 25% (12.1 

73 hours) of the data for this analysis was taken with the original configuration and 75% (44.1 hours) 

74 were taken with the new configuration. The observations were made in "wobble" mode in which 

75 the source is offset from the center of the field-of-vie w of the cameras to maximize efficiency in 

76 obtaining both source and background measurements (IFomin et al.lll994r). All dat a were processed 

77 and analyzed with the VEritas Gamma-ray Analysis Suite (VEGAS. ICoganl (12007)) wit h all results 
7B being successfully cross-checked by a second independent analysis package (lDanielll2008l ) . The event 
79 selection criteria (cuts) used are based on the image morphology (Mean Scaled Width and Length), 
so and the angular distance between the reconstructed position of the source in the camera plane 
si and the a priori known source location (lAharonian et al.lll997l ) and were optimized on a simulated 

82 source with a soft energy spectrum (r=4) and a 7% Crab Nebula flux (at 200 GeV). 

83 From the entire observational sample, a total of 822 events was de tected in exce ss of the 

84 estimated background, corresponding to a statistical significance of 6.4a (|Li &: Mal ll983). A two 

85 dimensional Gaussian fit to the VERITAS excess (VER J0416+011) is consistent with a point 
source located at 04 h 16 m 54 s ± 51„ f ± 61,., +1°06'35" ± l'47" tat ± 1'30" (J2000), and with the 



optical AGN position (jHewitt & Burbidgdll993l ). VER J0416+011 is shown in Figure 1 along with 
the Fermi-LAT and optical positions. 



The differential photon spectrum (see Figure 2) obtained from the VERITAS observations is 



-3.4±0.5 stat ±0.3 S!/Scm -2 s - 



L TeV" 



90 well fit by a power law of the form (1.6±0.3 s t a t±0.8 sys ) x 10 x ( Q 3TeV ) 

91 with a reduced chi-squared value (x 2 /d.o.f) of 0.57 for 2 degrees of freedom. The VERITAS spectral 

92 parameters are consistent with those measured by HESS. 

93 Due to its relatively low transit path on the sky for VERITAS, the data on 1ES 0414+009 were 

94 taken over a small range of zenith angles from 30° to 40°, resulting in an energy threshold 1 of 200 

95 GeV. The integral flux, obtained from 1ES 0414+009 during VERITAS observations (derived from 

96 the spectral fit shown in Figure 2) was (5.2±l.l s t a £±2.6 sys ) x 10~ 12 photons cm -2 , or 2% of the Crab 

97 Nebula flux above 200 GeV. The construction of a lightcurve, binned on the timescale of VERITAS 

98 "dark runs" (intervals between full moons when skies are dark) shows no evidence for variability 

99 with a straight line fit yielding a reduced \ 2 of 0.55 for a chance probability of 88% (see Figure 3). 

100 The HESS observations, which indicated a 0.6% Crab Nebula flux, were made between October 

101 2005 and September 2009, overlapping only slightly with the VERITAS observations. While the 



2 The energy threshold is defined as the energy corresponding to the maximum of the product function of the 
observed spectrum with the collection area of the instrument. 
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102 VERITAS integral flux is lar ger than the HESS measu rement, when statistical and systematic errors 

103 are taken into account (see lAbramowski et al.l ((2012)), the two flux measurements are consistent 

104 with each other. However, given the current data we cannot rule out some long term variability in 

105 flux. 



108 



3. Multiwavelength Observations 
3.1. Fermi-LAT 



Fermi-LAT (lAtwood et al.ll2009l ) analysis was performed on all available photons from 0.3 to 

109 300 GeV accrued from the start of full Fermi-LAT science operations (August 4, 2008) until July 14, 

no 2011 (MJD 54682-55756), overlapping to a large degree with the VERITAS observations (January 

in 2008 - February 2011). The data were analyzed using the ScienceTools v9r23pl package available 

112 from the Fermi Science Support Center (FSScjE Standard data quality cuts were applied as rec- 

113 ommended by the FSSC, with only "diffuse" class photons being used for this analysis (events with 

114 a high probability of being correctly identified as photons) . An un-binned analysis was performed 
us on all suitable photons coming from a region of radius 15° centered around 1ES 0414 +009 with 

116 all 2FGL sources (sources identified within the Fermi-LAT 2 year catalog) lying within 15° of 1ES 

117 0414+009 being modeled as well. In addition, both Galactic and extragalactic contributions to the 
us diffuse background were modeled out using the P6_Vll_Diffuse response functions. The Fermi-LAT 

119 systematic uncertainties are conservatively estimated at 10% for the energy range under evaluation 

120 in this work. 

121 The source is detected in the 0.3-300 GeV energy range with a test-statistic (TS) of 97 cor- 

122 responding to a detection significance of approximately 9.8u. A lightcurve binned on a monthly 

123 timescale shows no evidence for variability in the Fermi-LAT energy range with a straight line 

124 fit yielding a reduced x 2 of 0.81 for a chance probability of 76%(see Figure 3). The source spec- 

125 trum is well fit by a power law of the form (6.3+ 1.5 s tat) xl0~ 13 ( 100( ^ MeV )~ 1 - 9±0 - lstat photons 

126 cm - 2 s MeV . For modeling purposes, spectral data points were obtained by dividing the data 

127 into five equally spaced bins in log(E) and requiring that each bin have a TS value of at least 4. 

128 These spectral points along with the associated fit to the Fermi-LAT data are shown in Figure 

129 2. As can be seen, the spectra from the two instruments (Fermi-LAT and VERITAS) connect 

130 to each other quite well, with the attenuating effect of the extragalactic background light on the 

131 observed TeV flux being readily apparent. It should be noted that the two spectra were taken over 

132 different (although overlapping), extended time periods and cannot be considered to constitute a 

133 simultaneous measurement of the GeV- TeV spectrum. 



3 http://fermi. gsfc.nasa.gov/ssc/ 

4 http://fermi. gsfc.nasa.gov/ssc/data/analysis/LAT_caveats. html/ 
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Using the ScienceTools function gtprobsrc, the excess seen in the Fermi-LAT data is consistent 
with a point source located at RA = 04 /l 16 m 51 s ± 8 s stat , D EC = +1°04 , 36" 



± 2*2" (J2000), 



consistent within errors with both the AGN position from iHewitt &: Burbidgd (|1993l ) and the 
VERITAS source location (see Figure 1). It should be noted that all Fermi-LAT analysis results 
presented here are consistent with those derived for the source 2FGL J0416. 8+0105 in the 2 year 
Fermi-LAT point source catalog. 



3.2. Swift-XBT 



141 
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Between 2006 and 2009, a total of 6.24 ks of Swift-XRT (jBurrows et al.1 |2005j) observations 
were taken on 1ES 0414+009 in "photon counting" (PC) mode. T he reduced data p roducts used 
for this study were downloaded from the UK Swift Science CenteiR (|Evans et alJliooa ). Calibrated 
source and background files, along wit h the appropr iate response files, were used as inputs to the 
spectral fitting package XSPEC vl2.6 (|Arnaudlll996l ). The spectrum between 0.3 and 10 keV was 
fit while allowing the neutral hydrogen (HI) column density to vary freely. The resulting HI column 
density was fit as (1 .7 ± 0.16) x 10 21 cm -2 , higher than the value of 0.85 x 10 



22 cm 2 quoted in 



Kaberla et al.l (|2005l ). The spectrum is well fit (C-statistics were applied with less than 5% of model 



realizations providing a better fit) by a simple, absorbed powe r law (accounting for the red shift 
of the target) using the Wisconsin photoelectric cross sections ([Morrison &: McCammonlll983l ) or 
the "zwabs" model in XSPEC. The best fit for this model is provided by a normalization of (6.1 
±0.3) x 10~ 3 ph cm _2 s — 1 at 1 keV and a photon index of -2.4+0.1. All of these results are consistent 
with an analysis of the individual pointings analyzed separately, presented in Abramowski (2012). 
The deabsorbed source flux is measured as (20.3 ±0.1) x 10 _12 erg cm _2 s _1 with no evidence for 
strong variability in the PC mode data presented here. 

It should be noted, however, that ~2200 seconds of Swift-XRT observations taken in windowed 
timing (WT) mode on MJD 55231 and 55241 in 2010 indicate a large increase in the 0.3-10 keV flux 
from 1ES 0414+009. This data is fit by a simple, absorbed power-law with a derived photon index 
entirely consistent with the index from the PC mode observations, while showing a 0.3-10 keV flux 
of 58.8 (±0.2) x 10 _12 erg cm -2 s , a factor of 2-3 higher than previous observations. Since the WT 
mode observations suggest a relatively short flaring episode, while the rest of the multi-wavelength 
data were taken over a much larger timescale, the WT mode observations are not included in the 
modeling effort in this work. However, we report these flare data as an indication of overall variable 
behavior in the system. 



5 http: / /www.swift. ac.uk/ 
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3.3. MDM 



166 1ES 0414+009 was observed at the 1.3m McGraw-Hill Observatory of the Michigan-Dartmouth- 

167 MIT (MDM) Observator^ on the southwest ridge of Kitt Peak, during the nights of January 26 - 30, 

168 2011. Standard V, R, and I filters were used. The data were bias-subtracted and flat-field corrected 

169 using standard ro utines in IRAF, and c omparative photometry was done based on comparison star 

170 magnitudes from (IFiorucci et al.lll998l ). For the construction of VRI SEDs, the magnitudes were 

171 corrected with the extinction coefficients Ay = 0.39, A# = 0.32, and Aj = 0.23, as provided by the 

172 NASA/IPAC Extragalactic Databas<|] ( NED). In addi t ion to the MDM observations , we include 

173 archival U,B,V,R, and I filter data from Fiorucci et al] (2004) and Ulmer et al. ( 1983 ): along with 

174 infrared J, H, K values taken from Chen et al. ( 20051 ). 



175 4. Mult i wave length SED and Modeling 

176 The broadband SED constructed from the observations of 1ES 0414+009 is shown in Figure 4. 

177 Modeling of this SED proceeded by applying three distinct models of emission (two purely leptonic 

178 and one lepto-hadronic hybrid). The resulting predictions of the gamma-ray emission from these 

179 models were then mo dified to account f or the attenuating effect of the extragalactic background 
iso light (EBL) using the iFinke et al] (fcoid ) EBL model. 



The purely leptonic models employed to fit the SED are based on the premise that the high 
energy (GeV-TeV) emission is produced by inverse-Compton scattering of low energy ambient 
photons by a population of relativistic jet electrons (whose primary synchrotron emission makes 
up the lower energy X-ray peak). The variations between the models are based upon the nature 
of the photon field which is up-scattered: in synchrotron self-Compton (SSC) models, the target 
photon field is made of the original primary synchrotron photons (from the lower energy peak); 
while in external Compton (EC) models the emission includes the SSC contribution as well as 
inverse-Compton emission from a target photon field external to the jet. 



m 



The SSC model e mployed here is that of iBottcher &; Chiang! ((2002J) described also in detail 



Acciari et al.l (|2009i ). In this model, the emission originates from a spherical distribution of 



191 relativistic electrons (of radius R) traveling along the jet axis with a Lorentz factor T, corresponding 

192 to a jet speed of /3rc. The jet viewing angle o j, s results in a Doppler boost factor given by D = 

193 (r[l — /?rcos© fc s ]) _1 . In order to reduce the number of free parameters, we choose 0ofc s = 

194 leading to T = D. 

195 In models such as these, particles are injected with some distribution in energy and are sub- 

196 sequently cooled (i.e. lose energy by radiative mechanisms) or escape the accelerating region and 



http://mdm.kpno.noao.edu/ 
5 http://ned. ipac.caltech.edu/ 
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197 no longer contribute to emission. The models under consideration here result in an equilibrium 

198 between these factors (i.e. acceleration and cooling/escape). Here we parameterize the front end 

199 of the model by injecting electrons into the spherical region with a power-law distribution given by 

200 Q(t) = Qo7 _9e bounded by the cutoff parameters "femimlemax- The escape time in the model is 

201 parameterized by rj which corresponds to an escape time t esc of the particles of t esc = r/R/c. The 

202 resulting equilibrium particle distribution in the emission volume corresponds to a kinetic power 

203 L e in relativistic electrons. The synchrotron emission is determined by a tangled magnetic field B, 

204 corresponding to a power L# in Poynting flux. Furthermore, the equipartition parameter e e B = 

205 L^/L e is evaluated during the modeling process to give an evaluation of how reasonable the model 

206 is physically; it is reasonable to expect that most models that fit HBL emission well would generate 

207 values for e e s within the range of ~0. 1-1.0. 

208 The EC model used in this analysis incorporates all the above SSC model attributes while 

209 adding an additional emission component due to inverse-Compton scattering off an external photon 

210 source (corresponding to a standing perturbation of the jet flow). Hence its modeling parameters 

211 are identical to the SSC parameters along with the addition of the blackbody temperature and 

212 energy density of the external photon field (Tec, uec)- 

213 In contrast to the two leptonic models above, the lepto-hadronic model used for this analysis 

214 assumes that the high energy emission is dominated by proton synchrotron and pion decay emission 

215 processes. An additional model component is generated from the electromagnetic cascades from 

216 7 — 7 absorption between multi-TeV photons (from charged-pion decay products) and low-energy 

217 synchrotron photons from ultrarelativistic leptons. Similarly to the SSC and EC models, the 

218 protons are injected with a power-law spectrum (parameterized by 72(7) = 7 _<? j> bracketed by the 

219 lower and higher energy cutoffs ^pmimlpmax) corresponding to a kinetic power L p in relativistic 

220 protons, which we compare with the power L e in relativistic primary leptons and the Poynting flux 

221 power Lb- All parameter values described here, as well as resulting equipartition values, are listed 

222 in Table 1. 

223 As can be seen in Figure 4, while both leptonic models fit the TeV spectrum reasonably well, 

224 neither provides a satisfactory fit to the observed broad GeV peak (in contrast to the relatively 

225 sharp synchrotron peak). The extremely wide separation between the synchrotron and gamma-ray 

226 peak implies a very large 7 of the electrons (see Table 1), which then implies a very low B-field 

227 in order to achieve the observed synchrotron peak. A large Doppler factor (D = 40) can then 

228 help reduce the overall power requirements in the electrons to maintain values close to reasonable 

229 equipartition conditions. 

230 We find that the lepto-hadronic model provides a more satisfactory fit to the observed SED. 

231 The primary proton synchrotron emission matches the GeV - TeV emission very well, while the 

232 lower end of the GeV spectrum is reproduced by emission from electromagnetic cascades. It is 

233 clear however, that this model still over/underpredicts some of the Fermi-LAT GeV points and is 

234 far from equipartition. It is possible that additional components may be required in this model to 
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235 more accurately describe the system. Additionally, the model predicts a relatively strong flux in 

236 the very hard X-ray/soft gamma-ray regime; an energy range which is completely unexamined by 

237 the observations described in this work. 



23B 5. Summary 

239 The VERITAS array of telescopes has detected the distant blazar 1ES0414+009 in 56.2 hours 

240 of observations, resulting in a detection at a statistical significance of 6.4o\ The differential photon 

241 spectrum from the source is well fit by a power law with normalization and photon index parameters 

242 given by (1.6±0.3 stai ± 0.8 sys ) x 1CT 11 x _^- 3 ' 4±0 - 5 "»' ±0 - 3 '»» cm -2 s -i TeV -i i xhe inte gral flux 

243 from the source is measured as approximately 2% of the Crab Nebula flux above 200 GeV. While 

244 this is a larger integral flux than the HESS measurement (0.6% Crab above 200 G eV), when 



245 statistical and systematic errors are taken into account (see lAbramowski et al.l (|2012i )). the two 

246 flux measurements are consistent with each other. 

247 We have also presented multiwavelength data on 1 ES 0414+009 from Fermi-LAT , Swift— XRT, 
24B and the MDM observatory. We combine these data with archival optical data to construct a 

249 broadband SED which shows a very narrow synchrotron peak in contrast to the relatively broad 

250 GeV-TeV gamma-ray peak. The SED is not well fit by homogenous leptonic one-zone models, 

251 however a hybrid lepto-hadronic model provides a reasonable fit to the observed data in a strongly 

252 magnetically dominated jet. 

253 Although the lepto-hadronic model provides a better fit to the observed SED, it still does not 

254 accurately predict the apparently complicated form of the GeV-TeV emission. It should be noted, 

255 however, that both the GeV and TeV spectral data points have relatively large error bars and do 

256 not take systematic errors into account. Therefore, it is difficult to accurately constrain or improve 

257 any SED models until much more precise measurements of the GeV- TeV regime emission are made. 

258 The study of sources at such a relatively high redshift can shed light on more fundamental sci- 

259 ence topics such as the characteristics of the extragalactic background light (EBL) and intergalactic 

260 magnetic fields (IGMF). While we do not make any statements about these topics in this current 

261 work, as the catalog of high redshift TeV detections increases (with an accompanying increase in 

262 statistics on individual sources) our understanding of the properties of the EBL and IGMF will 

263 increase as well. A subsequent study of this source with these goals in mi nd will be the su bject of 



264 a future publication. Additionally, the next generation IACT array CTA (jActis et al.ll201ll ) (along 

265 with additional data collection with Fermi-LAT to improve the MeV-GeV spectral measurement) 

266 will be able to much more precisely measure the high energy emission of distant objects such as 

267 1ES 0414+009, providing valuable insights into the nature of some of the most distant TeV targets 
26B known as well as fundamental measurements of the EBL and IGMF. 



269 
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Parameter 


ssc 


EC 


Lepto-Hadronic 


Temin 


2xl0 5 


9xl0 4 


4.8xl0 3 


"lemax 


5xl0 6 


3xl0 6 


1.0x10 s 


Qe 


3.5 


3.5 


3.6 


Ipmin 


- 


- 


l.OxlO 3 


^pmax 


- 


- 


l.OxlO 11 


Qp 






1.8 


V 


120 


40.0 


3 


B(G) 


0.008 


0.044 


30 


r 


40 


40 


20 


R (cm) 


2.1xl0 17 


7xl0 16 


l.OxlO 16 


e o6s (°) 


1.43 


1.43 


2.86 


Strain. variability {h>TS) 


62.5 


20.9 


5.95 


Tec(K) 




1000 




u E c(erg/cm~ 3 ) 




3xl0" 9 




L e (erg/s) 


3.07xl0 44 


6.93xl0 43 


1.91xl0 41 


L B (erg/s) 


1.69 xlO 43 


5.69xl0 43 


1.35xl0 47 


e 


0.055 


0.82 


7.06xl0 5 


Lp 






8.0 xlO 43 


€ p B=Lb/L p 






1.6xl0 3 


£pe — Lg/Lp 






2.4xl0" 3 



Table 1: Summary of the emission zone parameters used for modeling the synchrotron self-Compton, 
external Compton, and lepto-hadronic models as described in the text. 
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Right Ascension (Deg) 



Fig. 1. — The two-dimensional significance map of the 1ES 0414+009 region from VERITAS 
observations made at VHE gamma-ray energies with the color scale representing units of standard 
deviation of the correspon ding excess. The cross represents the optical position of the AGN from 



Hewitt fc Burbidgd (|1993l ). the solid circle represents the best fit VERITAS position with associated 



statistical and systematic errors, and the dotted circle represents the position of the excess observed 
by Fermi-LAT (statistical error only). The white circle in the lower left corner represents the scale 
of the VERITAS point spread function of 0.12°. 
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Fig. 2. — The derived spectral points from both Fermi-LAT and VERITAS observations (non- 
simultaneous) of 1ES 0414+009. The dotted line shows the power law fit to the Fermi-LAT spec- 
trum extrapolated to TeV energies, while the solid line shows the power-law fit to the observed 
VERITAS TeV spectrum. 
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Fig. 3. — The VERITAS (top) and Fermi-LAT light curves for the observations detailed in this 
work. Errors shown on both light curves are statistical only. 
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Fig. 4. — The broadband SED constructed from the data sets described in the text. The top 
panel shows the SSC model results with contributions from primary synchrotron (dashed line), 
synchrotron self-Compton (dotted line), and EBL correct total (solid line). The middle panel 
shows the external Compton model results with contributions from primary synchrotron (dashed 
line), synchrotron self-Compton (dotted line), external Compton (dotted-dashed line), and EBL 
corrected total (solid line) . The bottom panels shows the lepto-hadronic results with contributions 
from electron synchrotron (dotted line), proton synchrotron (dotted-dashed), proton synchrotron 
with electromagnetic cascades (dashed line), and EBL corrected total. 



